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ABSTRACT 

We present observations of the unusual microlensing event OGLE 2003-BLG- 
235/MOA 2003-BLG-53. In this event a short duration (~7 days) low amplitude 
deviation in the light curve due a single lens profile was observed in both the MOA 
and OGLE survey observations. We find that the observed features of the light curve 
can only be reproduced using a binary microlensing model with an extreme (planetary) 
mass ratio of 0.0039+07 for the lensing system. If the lens system comprises a main 
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sequence primary, we infer that the secondary is a planet of about 1 .5 Jupiter masses 
with an orbital radius of ~ 3 AU. 

Subject headings: Gravitational lensing - (stars:) planetary systems - Stars:individual 
(OGLE 2003-BLG-235, MOA 2003-BLG-53) 



1. Introduction 

Gravitational microlensing occurs when a foreground object passes through or very near the 
line of sight of a background source star generating a well known symmetric light curve profile. 
If the foreground lens object is a star with an orbiting planet, then the presence of the planet may 
be detectable via a brief disturbance in the single lens light curve (Mao & Paczyhski 1991 ; Gould 
& Loeb 1992). This effect can potentially be utilized to detect planets with masses ranging from 
those of gas giants right down to terrestrial planets (Bennett & Rhie 1996). 

The short timescales of these deviations, ranging from a few days for giant planets to hours for 
terrestrial planets, and their unpredictability, present considerable challenges in any observational 
program. While some encouraging results have been obtained (Bennett et al. 1999; Albrow et al. 
2000; Rhie et al. 2000; Bond et al. 2002; Jaroszyhski & Paczyhski 2002), no firm detections of 
planets by microlensing have previously been obtained. 

In this Letter we report observations, obtained by OGLE and MOA, of the event OGLE 2003- 
BLG-235/MOA 2003-BLG-53 (hereafter 0235/M53) that was independently detected in both 
survey programs. We observed a 7 day deviation that was strongly detected in both surveys. We 
show that an extreme mass ratio binary microlensing model best reproduces the observed features 
in the light curve. 



2. Observations 

Presently, Galactic bulge microlensing events are discovered and then alerted by the two in- 
dependently operating survey groups OGLE (Udalski 2003) and MOA (Bond et al. 2001). The 
microlensing event OGLE 2003-BLG-235 (a=18:05: 16.35, £=-28:53:42.0, J2000.0) was first 
identified and alerted by the OGLE EWS system (Udalski 2003) on 2003 June 22. It was inde- 
pendently detected by MOA on 2003 July 21 and alerted as MOA 2003-BLG-53. 

OGLE observations were carried out with the 1.3 m Warsaw telescope at Las Campanas Ob- 
servatory, Chile, which is operated by the Carnegie Institute of Washington, equipped with a mo- 
saic CCD camera with 8192x8192 pixels. The images were obtained in the I band with an expo- 
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sure time of 120 seconds each. The observations presented here come from the OGLE-III phase of 
the OGLE survey and started in August 2001 . Additional photometry of the star was also collected 
during OGLE- II phase (1997-2000). This dataset indicates, however, no variability of the object 
during that period and was not used in further analysis. 

MOA observations were carried out from the Mt John Observatory in New Zealand with a 0.6 
m telescope equipped with a mosaic CCD camera with 4096x6144 pixels. The MOA images were 
obtained using 180 second exposures with a broad band red filter with its throughput centered on 
the standard I band. 

The photometry was derived using difference imaging analysis carried out independently by 
the OGLE and MOA teams on their respective datasets. This method is the current state-of-the-art 
for photometric accuracy in crowded fields (Alard 1999; Alcock et al. 1999). Our analysis resulted 
in two sets of time-series photometry in the I band corresponding to 183 OGLE measurements and 
1092 MOA measurements during 2000-2003. 

In Fig. 1 we show the light curve for this event on various timescales from March 2000 to 
the present. The long term behavior is typical for single point mass microlensing events, and it is 
similar to the almost 2000 other events discovered in the Galactic Bulge since 1993. The unique 
feature of 0235/M53 is, however, a short duration deviation from the profile expected for a single 
lens, seen clearly in both datasets during 2003 July 14-21. Moreover, a spike, characteristic of 
those binary microlensing events where the source enters or exits a bounded "caustic" region in 
the magnification map projected on the sky, was observed and well covered by MOA on 2003 July 
21. This caustic region was crossed in 12% of the overall lens Einstein radius crossing time. This 
short duration, combined with the small (~25%) amplitude of the photometric deviation in the 
caustic region interior, suggests an extreme mass ratio binary system. 

As well as regular monitoring in the I band, several V band observations were obtained by 
OGLE at various magnifications of the event. These were not used in the microlensing modeling, 
but they were used to constrain the source and lens star properties. By plotting the linearized fluxes 
in the I and V bands against each other, a model independent measurement of the color index of 
the source star was determined. We obtained V — I = 1.58 ± 0.02. Using E(V — I) = 0.82 
mag for the interstellar reddening towards the source (Sumi 2003), the corrected color index of 
(V — I) = 0.76 ± 0.02 indicates a G type source star. 

3. Light curve modeling 

The modeling of the observed light curve of 0235/M53 was performed independently by three 
groups using different methods to generate numerical binary microlensing light curves (Bennett & 
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Rhie 1996; Mao & Loeb 2001; Rattenbury et al 2002), and all three found the solution that is 
presented in Fig. 1. The observable quantities for all microlensing events are the Einstein radius 
crossing time, t E , the impact parameter u (in Einstein radius units) of the source star trajectory 
with respect to the lens center-of-mass (c.o.m.), and the time, t , of the closest approach to the 
c.o.m. For binary microlensing events, one also measures the mass ratio, q = Mi/M 2 , the trans- 
verse separation, a, of the lens components, and the position angle, 0, of the binary with respect 
to the source-lens transverse velocity. For caustic crossing events, one also measures the ratio, 
p = 6*/9 E , of the apparent angular radius of the source star to that of the Einstein ring. In addition 
to these 7 physical parameters, there are two linear scaling parameters between the magnification 
and the flux units for each passband, giving a total of 1 1 parameters for the modeling. In our 
modeling procedure, we searched for local x 2 minima using minimization procedures that allowed 
all 1 1 parameters to vary simultaneously. Our light curve modeling also employed a surface limb 
darkening profile appropriate for a G type star assuming a metalicity that is approximately solar. 

In Table 1, we list the physical microlensing parameters for the best fitting model shown in 
Fig. 1. This model has x 2 — 1390.49 for 1267 degrees of freedom and an extreme mass ratio of q — 
0.0039, which is a strong indication that the secondary may be a planet. Since microlensing light 
curves generally allow a much more accurate determination of the secondary:primary mass ratio, 
than the absolute mass of either body, the most sensible way to distinguish planetary microlensing 
events from those due to binary systems is through a criterion based on the mass ratio parameter 
q. There is a well known minimum in the distribution of mass ratios for binary stars and planetary 
systems, which is known as the "brown dwarf desert'.' There are few systems known with 0.01 < 
q < 0.1 (Halbwachs et al. 2000; Chabrier 2003; Mazeh & Zucker 2002). Thus it is sensible to 
define the boundary between stellar binary and planetary binary microlensing events in between 
these q values to minimize any possible ambiguity. This leads to a criterion of q < 0.03 for a 
planetary microlensing event, and so 0235/M53 is clearly in the planetary event category. 

We have carried out a systematic search in parameter space to try to find sets of model pa- 
rameters that might explain the observed light curve with a larger mass ratio. Binary microlensing 
models with q > 0.1 that traverse a caustic curve in ~ 7 days have much larger magnifications 
inside the caustic curve than is observed for 0235/M53. These binary lens events also have much 
larger deviations from a single lens light curve before and after the caustic crossings. As the mass 
ratio is decreased, the best fit light curves approach the observed light curves, with much weaker 
caustic crossing deviations. In Fig. 2, we show a close-up of the 7 day deviation with the best fit 
planetary model, compared with the best non-planetary model with q > 0.03 and best fit single 
lens model. The non-planetary binary models and single lens models are strongly disfavored with 
fit x 2 values that are larger by A% 2 = 210.96 and Ax 2 = 650.96 respectively. In both cases, the 
X 2 improvement for the best fit model is quite significant in both the MOA and OGLE data sets 
(see Table 1). The failure of the non-planetary binary model can be seen in the Fig. 2 inserts. This 
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model predicts both stronger caustic signals and significant deviations 13-20 days after the second 
caustic crossing that are not consistent with the observations. This model also shows some discrep- 
ancies at magnifications < 2, but these are not as strongly excluded due to the higher photometric 
uncertainties at lower magnification. 

Also shown in Fig. 2 is a planetary model with an earlier caustic crossing and a larger planet 
mass ratio: q = 0.0070. This fit represents a distinct local minimum of the x 2 surface, and is 
disfavored by only A% 2 = 7.37 or ~ 2.7<r. This is not accounted for by our 1-er uncertainty on 
q = 0.0039. Therefore, we have increased the upper error estimate on the planetary mass ratio to 
0.001 1, so that the actual uncertainty in q will be bounded by our error estimates at the 3-a level. 
In Table 1 we also list the parameters for these alternative models. 

Finally, the first OGLE observation after the second caustic crossing indicates a magnification 
below all of the binary models shown by about 3.6a. Such an outlier is not unusual because the 
real photometric error distributions for crowded field photometry generally have larger wings than 
a Gaussian distribution. Three of the 183 OGLE measurements are outliers from the best fit by 
> 3(7, and 16 are outliers by > 2a. These outlier points do not appear to cluster in the vicinity of 
the planetary deviation. If this single data point did indicate a real light curve deviation, it could 
be explained by a small variation to the planetary microlensing model, such as a moon orbiting the 
planet, but there is no non-planetary model that could help to explain it. 

4. Further constraints on the source and lens 

Most microlensing events have only a single measureable parameter, t E , that constrains the 
lens mass, distance, and tranverse velocity with respect to the line-of-sight to the source. How- 
ever, time resolved observations of binary event caustic crossings resolve finite source star effects 
and partly remove these degeneracies (Alcock et al. 2000; Witt & Mao 1994; Nemiroff & Wickra- 
masinghe 1994; Gould 1994), by allowing a measurement of the Einstein angular radius given by 
6>| = (4GM lcns /c 2 )(L> source — Acns)/(-D S ourcc Acc- 
using the flux parameters of the microlensing fit, we obtained I = 19.70 ± 0.15 for the source 
star and I = 20.7 ± 0.4 for the blended component. This source star magnitude, plus the V-I 
color from Section 2, can be compared to the bulge color magnitude diagram of Holtzman et al. 
(1998), and this indicates that the source is probably a bulge star near the main sequence turn-off. 
To determine the angular radius of the source star we used the color-color relations of Bessell & 
Brett (1988) together with empirical relations between V — K and surface brightness derived from 
interferometry observations of nearby main sequence stars (van Belle 1999; di Benedetto 1998). 

We find 9* = 0.50 ± 0.05 p:as, which combined with our measurement of p, yields 9 E = 
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520 ± 80 //as. This yields the following relation between the lens mass and distance 
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where x = D\ ena / D somce . If we combine this relation with the mass luminosity relations of Kroupa 
& Tout (1997) for main sequence stars, and require that the lens luminosity at a given distance does 
not exceed the blend flux, we obtain an upper limit (90% confidence) of Aens < 5.4 kpc. Thus, if 
the lens is a main sequence star, it must be in the Galactic disk. 

In Fig. 3, we show Eqn. 1 together with the results of a maximum likelihood analysis based on 
our measurements of the Einstein ring and its characteristic crossing time. The likelihood function 
was calculated using the Galactic disk models of Han & Gould (1996). We then obtain with 90% 
confidence: Acns = 5.2+2 g kpc from which we infer the lensing system to comprise an M2-M7 
dwarf star of mass 0.36!ta28 M© with a giant planetary companion of 1.5+?;2 Mj (Jupiter masses). 
The planet is in a wide orbit with a transverse separation of 3. 01$;? AU. 

Another possibility for the lens is that it could be a remnant object such as a white dwarf, 
neutron star, or black hole. If the lens is a white dwarf with mass 0.6 M Q , Eqn. 1 would place it a 
distance of 5.5 kpc. In this case, the microlensing parameters would imply a 2.5 Mj planet orbiting 
the white dwarf with a transverse separation of 2.8 AU. 



In Section 3 we concluded that the observed light curve of 0235/M53 is best described by a 
binary lensing model with an planetary mass ratio of q = 0.0039. Our definition of the planetary 
nature of the secondary lens by means of the mass ratio is optimal when the mass ratio can be 
measured, but it is useful to consider other possible definitions. Another potential dividing line be- 
tween planets and brown dwarfs is the Solar metalicity threshold for sustained Deuterium burning 
at 13.6 Mj, although Deuterium burning itself has little relevance for planet formation. 

The situation in the case of 0235/M53 was helped by the measurement of finite source effects. 
If the lens is a main sequence star, then as shown in the previous section, it must be an M dwarf 
with a ~ 1.5 Mj planetary companion 1 . There is a non-negligible chance that the lens is a white 
dwarf and a much smaller chance that it is a neutron star, but in both cases, a planetary companion 
below the nominal 13.6 Mj threshold is required. Only in the unlikely case of a massive black hole 
primary, could the secondary be outside the range traditionally associated with a planet. 



The only other M dwarf star known to have planetary companions is Gliese 876 (Marcy et al 2001). 



5. Discussion 
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There are some prospects for follow-up observations of this event. Our measurements of 
the finite source effects imply a proper motion of the lens with respect to the source of /x = 
#e Ae = 3.1 ± 0.4 mas/yr. High resolution imaging carried out ~ 10 years from now with JWST 
or adaptive optics systems should be able to resolve the lens and source stars providing direct 
measurements (Han & Chang 2003; Alcock et al. 2001) of the color and brightness of the lens, as 
well as confirmation of the proper motion measurement. 

We present these observations as a demonstration of the planetary microlensing phenomenon. 
The power of microlensing is in its ability to acquire statistics on many systems (Bennett & Rhie 
2002). These include planets in wide orbits, very low mass planets, and even planets in other 
galaxies (Covone et al. 2000; Bond et al. 2002). The challenge now to the microlensing community 
is to develop effective strategies to find more planetary microlensing events. 

Numerical photometry of OGLE 2003-BLG-235/MOA 2003-BLG-53 is available from the 
websites for OGLE http://ogle.astrouw.edu.pl and MOA http://www.physics.auckland.ac.nz/moa. 
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B. Paczyhski. A.U., I.S., and K.Z. also acknowledge support from the grant "Subsydium Profe- 
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2820 2840 2860 2880 
HJD - 2450000 

Fig. 1 . — Light curve with best fitting and single lens models of 0235/M53. The OGLE and MOA 
measurements are shown as red filled circles and open blue circles, respectively. The top panel 
presents the complete dataset during 2003 (main panel) and 2001-2003 OGLE data (inset). For 
clarity, the error bars were not plotted, but the median errors in the OGLE and MOA points legend 
are indicated. The lower panel is as the top panel but with the MOA data grouped in 1 day bins 
except for the caustic crossing nights, and with the inset showing MOA photometry during 2000- 
2003. The binary and single lens fits are plotted in flux units normalized to the unlensed source 
star brightness of the best planetary fit, and they are indicated by the solid black and cyan dashed 
curves, respectively. 
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Fig. 2. — Light curve and modeling of 0235/M53 during the caustic traverse. The models shown 
are the best fits for various regimes of the lens system parameter space, normalized as in Fig. 1. 
These models are: single lens case (cyan, long dash curve), binary lens with q > 0.03 (magenta, 
short dash line), planetary lens with caustic entry before day 2835 (green, dotted line), and the best 
overall fit with q = 0.0039 (dark, solid line). The insets show the second caustic crossing and a 
region of the declining part of the light curve where the best fit non-planetary binary lens model 
fails to fit the data. MOA data on days other than the caustic entry and exit (days 2835 ±0.5 and 
2842 ±0.5) are placed in one day bins. 
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Fig. 3. — Constraints on the distance to 0235/M53 and its mass. Our measurement of the angular 
size of the Einstein ring constrains the lens mass and distance to lie on the solid curve. The curve is 
bounded above and below by the dark dashed curves, due to the uncertainty in this measurement. 
The likelihood function (shown by the magenta curve) is for a main sequence disk star prior for 
the lens star with the sharp cut-off being due to the constraints imposed by the measured lens flux. 



Table 1 . Microlens Model Parameters 



Model 


M p /M* 




a proj/^?E 


<t> 


u 


to 


*E 
(days) 


^source 

mag 


x 2 

(1267 dof) 


Xmoa 
(1089 dof) 


XOGLE 

(178 dof) 


Best Fit 


0.0039(1^) 


0.00096(11) 


1.120(7) 


223?8(1?4) 


0.133(3) 


2848.06(13) 


61.5(1.8) 


19.70(15) 


1390.49 


1151.00 


239.50 


Early Caustic 


0.0070 


0.00104 


1.121 


218?9 


0.140 


2847.90 


58.5 


19.62 


1397.87 


1149.37 


248.49 


Best Non-planet 


0.0300 


0.00088 


1.090 


187?9 


0.144 


2846.20 


57.5 


19.68 


1601.44 


1229.47 


371.98 


Single Lens 










0.222 


2847.77 


45.2 


19.10 


2041.45 


1624.17 


417.28 



Note. — 



The units for to are HJD - 2450000. 



